Stochastic Resonance Phenomena in Quantum Many-Body Systems 
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We discuss stochastic resonance (SR) effects in driven coupled quantum systems. We construct 
dynamical and information theoretic measures of the system's response that exhibit a non-monotonic 
behaviour as a function of the noise strength. We analyze the relation between lack of monotonicity 
in the response and the presence of quantum correlations and identify parameter regimes where the 
breakdown of a linear response can be linked to the presence of entanglement. We also show that 
a chain of coupled spin systems can exhibit an array-enhanced response, where the sensitivity of 
a single resonator is enhanced as a result of nearest-neighbour coupling. These results enlarge the 
domain where SR effects exist and should be observable in arrays of superconducting qubits. 

PACS numbers: PACS numbers: 03.67.-a, 03.67.Hk 



Recently, there has been an increasing interest in the 
observation of entanglement and correlation phenomena 
in arrays of super-conducting and other solid state qubits 
realizations In general, the presence of noise and 
finite temperatures is considered detrimental and noise 
levels are thus aimed to be minimized. Here, however, we 
will explore situations where it is advantageous to main- 
tain a finite, not necessarily minimal, noise strength. 
The response of an open quantum system to a weak pe- 
riodic forcing can exhibit a resonance-like dependence on 
the noise strength A canonical example to illustrate 
this form of quantum SR ( QSR ) is provided by the pe- 
riodically driven biased spin-boson model H. Despite 
quantum coherence was believed to contribute to the dis- 
appearance of SR effects [3| , recent work has shown that 
QSR should also be displayed by systems whose dissi- 
pative dynamics obeys conventional Bloch equations [B} 
and different experimental realizations in quantum op- 
tics have been proposed @. We will show that noise- 
enhanced effects are also present in the steady state re- 
sponse of quantum-mechanically correlated systems and 
quantify the presence of stochastic resonance in terms of 
both dynamical and information theoretic measures. 
Our system consists of an array of N driven, longitu- 
dinally coupled spin-1/2 systems (qubits). The system 
is subject to a noisy environment modelled by a set of 
harmonic oscillators such that each qubit couples trans- 
versely to its own bath. The global Hamiltonian is given 
by, 
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where % — 1, X 1 = J2k^ k ( a k + a k) denotes the bath's 
force operator and a\ = |1),(0|. The external driving 



is parameterized by its intensity, as given by the Rabi 
frequency f2j, and the detuning from the qubit frequency 
Si = uj'q — io L [7| . We will consider situations where the 
driving is weak and the external Rabi frequency is smaller 
than the interqubit coupling, Q < J [8j. Within the 
rotating wave approximation and for a Markovian bath, 
we obtain an effective Hamiltonian for the ./V-qubit array, 
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where, H coh = - YhLi % a \ 

f^j cf 1 x is the coherent part of the Hamiltonian in 
the interaction picture. This master equation treatment 
is valid in the parameter regime Qi/u> <C l,Tifi/u> <C 
l,8i/u> <C 1 and J/uj <C 1, where u — mining, uj c } for a 
suitable frequency cut off lo c of the bath, and all expres- 
sion in this work will then be correct to lowest non-trivial 
order in fl; , T.; , Si and J . 

The noise strength on qubit i at a temperature T is given 
by the product Ti n, where the explicit functional form 
of the decay rate Ti depends on the spectral properties 
of the bath and n denotes an effective boson number 
that depends on the bath's temperature T; both pa- 
rameters are, in principle, controllable. For instance, 
for superconducting qubits, the parameter T may be a 
measure of the fluctuations in the gate voltages which, 
if desired, may be amplified on demand [10]. As it will 
be clear below, our argument is that the steady state 
response of the system, as quantified by different fig- 
ures of merit, will be optimized at intermediate noise 
levels and therefore, trying to reduce the environmental 
noise to as small as possible values, does not necessar- 
ily provide an optimal universal strategy to maximize 
coherent effects. Let us consider first the case where 
N = 2 and T = 0. By integrating the time evolution 
f> = -i[H cS ,p] + J2iFt(n + \)o l _pa\ + J2i <?+pcr!_, 
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FIG. 1: Entanglement time evolution for two weakly driven 
qubits with longitudinal coupling of strength J/Q — 1.5 at 
T = 0. If the noise strength F is sufficiently large, the system 
is inseparable in the steady state (pink dashed line). 



for fixed values of the coupling J and the driving Q = 
fl i = Q2, we can analyze whether our system of weakly 
driven qubits, initially prepared in their ground state, 
develops quantum correlations in time. We employ the 
entanglement of formation Ep to quantify bipartite en- 
tanglement [ll[ . In Figure 1 we observe that the system 
will be entangled in the steady state only for certain fi- 
nite values of T. Perhaps surprisingly, it is the largest 
value of the noise strength the one that yields steady- 
state entanglement. The steady state of the system can 
be computed analytically. Calling r = T/il, s — J /ft and 
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where k = 3 + 2r 2 + t 2 + Ar 2 s 2 and . . . refer to the suit- 
able complex conjugate matrix element. The system is 
entangled, and have a negative partial transpose, only if 
r > T t h, where 
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is the noise threshold. If T < T t h, the state is separa- 
ble. This behaviour is illustrated in Figure 2 where the 
dashed line corresponds to the bipartite entanglement in 
the steady state as quantified by the Ep as a function 
of the noise strength T. As a result of the constraint 
given by Eq.([2]), any entanglement measure exhibits an 
initial domain of vanishing entanglement for weak noise 
where the state is separable. The smaller the qubit in- 
teraction strength J, the larger the value for the noise 
required for the driven spins to be entangled. When T 
rises above threshold, the steady state entanglement in- 
creases monotonically up to a maximum at certain opti- 
mal noise strength and decreases steadily for higher val- 
ues of r 12j . This functional form for the bipartite en- 
tanglement in the domain where T > T t h is reminiscent of 



FIG. 2: Stochastic resonance phenomena quantified in terms 
of information theoretic measures for a systems of two coupled 
and weakly driven spins (J/Q = 1.5) at a zero temperature. 
Presented results are for 8 = but deviations up to 5/Q ~ 
10 -2 yield very small deviations from the exact resonance 
behaviour. 



stochastic resonance Q and it had been observed before 
in the context of incoherently driven quantum systems 
[l3j |. However, to argue for the system to display SR 
in the conventional sense, we need to construct suitable 
measures of the system's information content and show 
the characteristic non-monotonic behaviour as a function 
of the noise strength that is typical of an SR response. A 
possible information-theoretic measure is provided by the 
system's mutual information I12 = S\ + $2— S12, where S 
denotes the von Neumann entropy, S(p) = —tr(p log 2 p) 
[pi ] . The blue solid line in Figure 2 is the mutual informa- 
tion of the steady state p\ s 2 - Correlations increase mono- 
tonically up to a maximum corresponding to a certain 
optimal noise strength above which I12 decreases. This 
characteristic response is also obtained for the difference 
between the mutual information and the entanglement 
(pink solid line). Longer chains, as detailed later, also ex- 
hibit this type of non-monotonic response and we there- 
fore argue that SR can be observed beyond the purely 
incoherent regime analyzed in 
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Alternatively, stochastic resonance can be character- 
ized using a dynamical measure of the system's response 
to the external driving. In the case of single quantum sys- 
tems, it has been proposed to use as an appropriate figure 
of merit the expectation value of a suitable Pauli oper- 
15j . This suggests considering local observables of 
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the form Sn = ((1/N) Yli=i a l) (C = x i z ) to character- 
ize the system's dynamical response in the multipartite 
scenario. Note that Sn is easily accessible experimen- 
tally and may be even measured without local control. 
We can evaluate analytically the form of the signal £2 = 
((ovj, + a 2 )/2) = 4sr 2 /k, which is non-monotonic as a 
function of both the coupling J between qubits and the 
coupling to the reservoir, reaching a maximum at cer- 
tain intermediate values. The signal £2 is maximum for 
T = v2^) independently of J, with <S2, roa x = s /(2 + s 2 )- 
The amplitude of the signal is maximal for s — J/Q = \[2 
and as the ratio s increases, the response becomes weaker 



3 





I 1g E0F lg 


It \ 

If '»,»» 

// *'*■«*.--.- 





5 10 15 20 25 30 35 40 

Tx 0.05 (Q) 



FIG. 3: System's response for a chain of N = 4 qubits at 
T = with Ji/Qi — 1.5, (i = 1 ... 4) as quantified by an infor- 
mation theoretic measure, the quantum mutual information. 
Under the evolution given by -ffeff, entanglement is restricted 
to nearest neighbours only and Ef 13 = Ef 14 — 0. 



dynamical measures of SR are reduced with increasing T 
(See Figure 5 as an illustration for typical superconduct- 
ing qubits operating temperatures). The mutual infor- 
mation Iij decreases monotonically with n, as illustrated 
by the blue solid line in Figure 4 for qubits 1 and 2 in a 
chain with N = 4, and so does the steady state entan- 
glement, which becomes zero for T sufficiently high. We 
should stress that I 12 quantifies 'total' correlations, not 
distinguishing between quantum and classical contribu- 
tions. Given that Iij measures total correlations, and the 
Ep is a measure of entanglement (but not the only one), 
their difference ly — Epij may be viewed as a reasonable 
attempt to measure the classical correlation content of a 
quantum state. This quantity displays the characteristic 
SR shape as T increases, with the maximum response 
being reached at the point where Ep 12 = (pink solid 
line in Figure 4) (ljj . 



and loses its resonance-like shape, resembling the shape 
of a saturation curve . In general, for arbitrary values of 
f2 and J, the transition from linear to non-linear response 
has no direct relation with the presence of quantum cor- 
relations in the system. The system will be separable in 
the linear region if T t h > V2D, while it can be entangled 
and responding linearly if Tth < T < \/2il. However, 
tuning the nearest neighbour coupling J to the value 
yields a maximal response at exactly the transi- 
tion point from PPT i— ► NPT states, i.e. from separable 
to entangled subsystems. As a result, if ft < \/8J, we 
can ensure that the breakdown of the monotonic response 
when varying T implies the presence of entanglement in 
the system. This happens in the regime where Q < J 
for both the dynamical response £2 and the SR infor- 
mation theoretic measures: The maximum values of the 
response are reached for a noise strength above thresh- 
old and the lack of monotonicity allows to conclude that 
there is entanglement in the system. Note, however, that 
the opposite conclusion, i.e. monotonicity ==> separabil- 
ity, is not correct in general, even for weak driving. 

SR phenomena, as quantified by information theoretic 
and dynamical measures, should also be observable in 
chains of longitudinally coupled weakly driven spin sys- 
tems. With the type of interaction we have considered, 
quantum correlations are confined to nearest neighbours 
and retain the same qualitative behaviour discussed in 
the N = 2 case, with quantum correlations only arising 
above a noise threshold that decreases as function of N. 
In Figure 3 we have characterized the system's response 
for a chain with N = 4 in terms of the quantum mutual 
information between qubits 1 and j, (j = 2, 3, 4). We ob- 
serve a resonance-like response as a function of the noise 
strength at T = for both nearest and distant neigh- 
bours. So far we have considered the environment to be 
at zero temperature. All phenomena described so far are 
robust in the presence of a finite T, with the net result 
that maximum values in both information theoretic and 





FIG. 4: Information the- FIG. 5: System's response for 

oretic measures as a func- a chain of N = 4 qubits with 

tion of the bath's mean bo- J/Q = 1.5 as quantified by 

son number (temperature) £4= (2 =1 a% ^l^) as a func- 

for a chain of N=4 spin with tion of the noise strength T 

Ji/Qi = 1.5 and Ti/Qi = 1, and for different values of the 

i = 1 . . . 4. external temperature. 

The non-monotonicity of the response is also apparent 
in the dynamical measure 1S4, plotted in Figure 5 for 
increasing values of the mean thermal boson number n. 
At a given T, the value of T that maximizes the response 
is now a function of both f2 and J and a numerical 
analysis shows that we can link entanglement and lack of 
monotonicity whenever outside the regime where £1 3> J. 
It is interesting to compare this global response with 
the individual signal that can be obtained from each 
qubit alone. Classical stochastic resonators are known 
to display an array-enhanced SR, where the collective 
dynamics yields the amplification of some suit ably 
defined signal-to-noise ratio of a single oscillator 17], [l8| • 
On exact resonance, the steady state expectation value 
(a x ) for a driven isolated qubit transversely coupled to 
a bosonic environment at zero temperature is strictly 
zero, given that (a x ) = poi + Pro ~ For a finite 

detuning, the single qubit response is a monotonically 
decreasing function of the noise strength. Specifically, 
I (a x ) ss |= flS/(S 2 + (r/2) 2 + (n/2) 2 ). Introducing 
interqubit coupling leads to a build up of a real part 
in the single qubit coherences which is proportional 
to the coupling strength J to first order. For N = 2, 
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FIG. 6: Single qubit steady state response for weakly driven 
arrays of 2, 4 and 5 qubits with J/Q. = 1.5 and n = 0.01. 
The zero signal corresponds to a resonantly driven uncoupled 
oscillator. Interqubit coupling yields an enhanced response as 
a function of the number N of oscillators in the array. 



(a x ) = 4 /rT a and the single qubit response displays 
the typical SR profile, while increasing the size of the 
array yields a sharper signal with a maximum value that 
increases as a function of N and that is obtained for 
increasingly smaller values of the noise parameter, as 
noted in Figure 6. The dashed-dotted line corresponds 
to the signal for the first qubit in an array of N = 4, 
while the dashed line corresponds to the response of the 
first qubit alone in an array of 5 spins. If we compare the 
numerical value of the single qubit expectation values 
with the global response specified in Figure 5, we note 
that the single qubit response is enhanced as a result 
of the coupling, which shows the persistence of array 
enhancement effects in chains of quantum spin systems. 
It seems clear from the present analysis that stochastic 
resonance phenomena predicted for coupled classical 
resonators, including those obeying an Ising model 
[3], may also be observed within coupled qubits and 
coexist with the presence quantum correlations. These 
results are amenable to experimental verification on a 
variety of proposed qubit realizations [2(J. Solid state 
architectures, and in particular superconducting qubits 
[Io| , combine noisy environments with a degree of qubit- 
environment tunability that made them particularly 
suited to demonstrate SR in the terms discussed in this 
Letter. 
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